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TROPICAL  CYCLONE  WIND  PROBABILITY  FORECASTING 


FOR  THE  EASTERN  NORTH  PACIFIC  (EPWINDP) 

1.0  Introduction 

The  concepts  and  development  of  tropical  cyclone  wind 
probability  forecasting  for  the  eastern  North  Pacific  follow 
closely  those  originally  presented  by  Jarrell^  (1981)  for 
the  western  North  Pacific.  This  is  the  third  in  a  series  of 
wind  probability  forecasting  program  reports  covering  differ¬ 
ent  ocean  basin  areas.  The  initial  report  on  the  western 
North  Pacific  was  followed  by  a  similar  report  on  the  western 
North  Atlantic  (Jarrell^,  1981).  This  report  and  the  two 
previous  are  based  upon  an  extension  of  the  concepts  developed 
for  tropical  cyclone  strike  probability  forecasting 
(Jarrell^,  1978).  As  these  concepts  have  been  previously 
presented  and  the  programs  are  presently  operational,  a  re¬ 
development  of  the  concepts  will  not  be  presented  in  this  re¬ 
port.  A  brief  summary  will  be  presented  as  a  background  for 
the  new  reader.  If  more  detailed  information  is  desired  the 
reader  is  referred  to  the  above  mentioned  documents.  Differ¬ 
ences  between  the  development  of  eastern  Pacific  tropical 

^Jarrell,  J.D.,  1981;  Tropical  Cyclone  Wind  Probability 
Forecasting  (WINDP),  NAVENVPREDRSCHFAC  Contractor  Report  CR  8- 
1-03. 

^Jarrell,  J.D.,  1981;  Atlantic  Hurricane  Wind  Proba¬ 
bility  Forecasting  (WINDP),  NAVENVPREDRSCHFAC  Contractor  Re¬ 
port  CR  81-04. 

^Jarrell,  J.D.,  1978;  Tropical  Cyclone  Strike  Probability 
Forecasting,  NAVENVPREDRSCHFAC  Contractor  Report  CR  78-01. 
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cyclone  wind  probabilities  and  that  of  the  Atlantic  is  the 
crux  of  this  report  and  these  will  be  described  in  detail. 

The  tropical  cyclone  wind  probability  forecasting  model 
generates  estimates  of  the  probability  of  30  and  50  knot  winds 
occurring  at  a  point  given  that  the  tropical  cyclone  occupies 
a  particular  position.  The  model  includes  many  features  of 
the  strike  probability  model  which  is  based  on  an  analysis  of 
position  forecast  errors  to  determine  the  probability  of  a 
tropical  cyclone  occupying  a  particular  geographic  position. 
Wind  profile  errors  in  the  forecast  of  maximum  wind  and  in  the 
forecast  radius  of  30  and  50  knot  winds  are  similarly  analyzed 
to  determine  the  30  and  50  knot  probabilities  of  occurrence. 

Jarrell^  (1978),  in  his  development  of  the  strike  prob¬ 
ability  model,  based  the  theory  of  strike  probability  on  three 
assumptions: 

1)  All  tropical  cyclone  forecasts  are  subject  to  error; 

2)  Difficulty  of  forecast  and  size  of  forecast  error  are 
related;  and 

3)  The  occurrence  of  errors  is  random  and  approximates  a 
multimodal  bivariate  normal  probability  distribution. 

Subsequent  studies  by  independent  investigators  in  three 
ocean  basins  frequented  by  tropical  cyclones  verified  these 


assumptions.  Nicklin^  (1977),  Thompson  and  Elsberry^ 
(1979),  and  Crutcher®  (1980)  using  data  from  the  western 
Pacific,  eastern  Pacific  and  the  Atlantic,  respectively  ,  also 
developed  individual  methods  to  group  forecasts  into  three 
classes  of  forecast  difficulty.  A  general  relative  classifi¬ 
cation  evolved  in  each  study  to  yield  forecast  error  groups  of 
below  average,  average,  and  above  average  errors.  These 
groups  were  also  referred  to  as  Class  I  (easy  forecasts). 
Class  II  (average  forecasts)  and  Class  III  (difficult  fore¬ 
casts)  in  previous  reports.  Each  investigator  also  utilized 
sufficient  statistical  data  and  prescribed  sufficient  parame¬ 
ters  to  describe  the  bivariate  normal  distributions  for  each 
of  the  three  classes  and  for  forecasts  of  24,  48  and  72 
hours. 


^Nicklin,  D.S.,  1977:  A  Statistical  Analysis  of  Western 
Pacific  Tropical  Cyclone  Forecast  Errors,  Naval  Postgraduate 
School,  M.S.  Thesis,  June. 

^Thompson,  W.J.  and  R.  L.  Elsberry,  1981:  An  Analysis  of 
Eastern  North  Pacific  Tropical  Cyclone  Forecast  Errors, 
Monthly  Weather  Review,  V  109,  pp.  1930-1938. 

®Crutcher,  H.L. ,  1980:  Tropical  Storm  Forecast  Error  and 
the  Bivariate  Normal  Distribution.  13th  Tech.  Conf.  on  Hurri¬ 
canes  and  Tropical  Meteorology,  AMS,  Miami,  FL. 


2.0  Model  Description 

2.1  Basis  for  Utilization  of  Atlantic  Model 

Analysis  of  eastern  Pacific  tropical  cyclones  has  been 
hampered  by  both  a  lack  of  historical  data  from  the  region  and 
a  lack  of  targets  (land  mass  or  ships)  which  would  provide  the 
necesary  incentive.  Therefore  the  wind  characteristics  of  the 
eastern  Pacific  cyclones  have  not  been  as  well  studied  as 
those  of  the  western  Pacific  and  the  Atlantic. 

A  study  by  Upton^  (1973)  indicated  that  eastern  Pacific 
tropical  cyclones  may  be  more  closely  related  to  Atlantic 
hurricanes  than  western  Pacific  storms  in  extent,  intensity 
and  other  features.  His  value  of  23.5  n  mi  for  the  radius  of 
maximum  wind  for  eastern  Pacific  cyclones  would  place  them 
well  within  the  15-30  n  mi  range  of  Shea  and  Gray®  (1972) 
for  moderate  Atlantic  storms.  Reports  by  Frank  and  Clark® 
(1979)  and  Simpson  et  al^®  (1968)  indicate  that  the 

^Upton,  T.G.,  1973:  An  Analysis  of  the  Thermal  and  Circu¬ 
lation  Features  of  Eastern  North  Pacific  Cyclones  Using  Air¬ 
craft  Reconnaissance  Data.  Naval  Postgraduate  School,  M.S. 
Thesis . 

®Shea,  D.J.  and  M.W.  Gray,  1972:  The  Structure  and 
Dynamics  of  the  Hurricane's  Inner  Core  Region.  Colorado  State 
Univ.,  Atmospheric  Science  Paper  182,  Ft.  Collins,  CO. 

®Frank,  N.L.  and  G.  Clark,  1979:  Atlantic  Tropical  Sys¬ 
tems  of  1978,  Monthly  Weather  Review,  V  107,  pp.  1035-1041. 

^®Simpson,  R.H.  et  al ,  1969:  Atlantic  Tropical  Distur¬ 
bances  of  1968,  Monthly  Weather  Review,  V.97,  pp.  240-255. 


eastern  North  Pacific  tropical  cyclones  may  be  considered  a 
westward  extension  of  the  North  Atlantic/Caribbean  tropical 
cyclone  area. 

More  importantly  perhaps  is  the  implication  by 
Jarrell (1981)  that  Classes  1,  2  and  3  group  similari¬ 
ties  (or  differences)  in  the  different  ocean  basins  may  be  due 
in  part  to  collection  methodologies  and  forecast  emphasis  and 
techniques.  Longer  range  (72  hr)  forecasts  are  much  better  in 
the  western  Pacific  than  in  the  Atlantic  while  short  range 
forecasts  are  better  in  the  Atlantic.  The  western  Pacific 
storms  are  forecast  by  the  USN/USAF  Joint  Typhoon  Warning 
Center  located  on  Guam  and  emphasis  is  on  long  range  forecasts 
related  to  military  oriented  operations  and  decision  making. 
Atlantic  hurricanes  and  eastern  Pacific  tropical  cyclones  are 
forecast  by  the  National  Weather  Service  whose  primary  in¬ 
terest  is  short  range  public  warning  forecasting.  This  would 
seem  to  provide  further  rationale  to  use  the  Atlantic  model 
for  the  eastern  Pacific. 


^^Jarrell,  J.D.,  1981;  Atlantic  Strike  Probability 

Program,  NAVENVPREDRSCHFAC  Contractor  Report  CR  81-04. 


2.2  The  Eastern  Pacific  Wind  Probability  Model 

Tropical  cyclone  wind  probability  was  developed  as  an  ex¬ 
tension  of  the  western  Pacific  strike  probability  program  by 
Jarrell^  (1981).  The  concepts  of  the  strike  portion  of  the 
eastern  Pacific  model  are  similar  to  those  developed  for  the 
western  Pacific.  However,  the  eastern  Pacific  model  uses  a 
discriminant  analysis  routine  developed  specifically  for  that 
region  (Thompson  and  Elsberry  )  and  the  Tsui-*-^  routine 
(used  in  the  Atlantic  wind  probability  model)  to  determine 
wind  profiles  and  to  handle  asymmetrical  storm  distribution. 

The  eastern  Pacific  strike  probability  program  used  the 
UCLA  BIOMED  discriminant  analysis  routine  (Dixon^^,  1975) 
to  develop  functions  from  eastern  Pacific  data  to  discriminate 
on  forecast  error.  Predictands  were  forecast  error  Classes  1, 
2  and  3.  In  contrast  to  western  Pacific  results,  there  was 
generally  poor  discrimination  between  Classes  1  and  2  (similar 
to  Atlantic  results)  while  Class  3  appeared  to  be  well  sepa¬ 
rated  from  1  and  2  (again  similar  to  the  Atlantic  results). 
Table  1  depicts  average  forecast  errors  for  classes  1,  2  and  3 
for  each  ocean  basin. 


^^Tsui,  T.L.,  1980;  Surface  Wind  Distribution  of 
Western  North  Pacific  Tropical  Cyclones;  13th  Tech.  Conf.  on 
Hurricanes  and  Tropical  Meteorology;  Miami,  FL. 

^^Dixon,  W.J. ,  1975;  BMDP  Biomedical  Computer  Programs, 
University  of  California  Press,  Berkeley. 
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TABLE  1 


Average  forecast  errors  <NM)  for  Class  1,  2  and  3  tropical 
cyclones  for  western  Pacific  (WP),  eastern  Pacific  (EP)  and 
North  Atlantic  (NA)  basins. 


WP  EP  NA 


1 

2 

3 

1 

2 

3 

1 

2 

3 

24 

hr 

99 

130 

148 

94 

97 

132 

81 

99 

160 

48 

hr 

204 

251 

286 

176 

188 

154 

200 

236 

325 

72 

hr 

324 

378 

407 

275 

297 

393 

362 

394 

477 

If  we  interprete  poor  discrimination  as  "nearly  equal  forecast 
errors  for  different  classes",  we  can  compare  the  classes  for 
respective  ocean  basins.  Averaging  differences  between 
Classes  1  and  2  (24,  48  and  72  hr)  and  then  between  Classes  2 
and  3  for  each  basin,  we  obtain  the  following: 


Differences: 

diff 

EP 

diff 

Class  1  to  Class 

2 

44.0 

31.7 

12.3 

16.3 

28.6 

Class  2  to  Class 

3 

27.3 

38.3 

65.6 

12.0 

77.6 

The  average  error  difference  (in  NM )  between  classes  can  be 
(depicted  as  diff)  reinforces  the  arguments  of  section  2.1; 
that  is,  the  eastern  Pacific  classes  seem  more  closely  related 
to  the  Atlantic  (either  in  nature  or  forecast  emphasis)  than 
the  western  Pacific. 
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Thus  while  using  the  unique  eastern  Pacific  strike  proba¬ 
bility  program,  this  model  uses  the  wind  probability  deriva¬ 
tion  and  concepts  used  in  the  Atlantic. 

The  Atlantic  wind  profiles  and  tropical  cyclone  asymmetry 
are  derived  using  methods  produced  by  Tsui^^.  Tsui,  using 
wind  radius  data  from  tropical  cyclone  warnings  over  a  12-year 
period  (1966  to  1977),  determined  that  the  profile  of  the  tan¬ 
gential  wind  speed  along  the  radial  axis  was  exponential.  He 
further  determined  that  maximum  wind  and  persistence  could  be 
statistically  related  to  the  size  of  the  storm  and  that  the 
asymmetric  shape  of  the  storm's  wind  pattern  could  be  corre¬ 
lated  to  the  forward  speed  of  movement  of  the  storm.  A  simple 
empirical  relationship  was  derived  to  provide  an  estimation  of 
any  wind  radii: 

V/Vmax  =  (-0.693R), 

where  V  is  the  wind  speed  of  interest,  V„  is  the  maxi- 

mum  wind  speed  of  the  storm,  and  R  is  a  ratio  of  the  radius 
associated  with  V  to  the  radius  associated  with  one-half  of 
^max  ^^half^  respectively.  Asymmetrical  storm  con¬ 

figuration  is  accommodated  in  the  profile  by  an  empirical 
adjustment  of  the  one-half  radius,  ^half  ’ 
bearing  relative  to  direction  of  motion  and  translation 
speed . 
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3.0  Testing  the  Eastern  North  Pacific  Wind  Probability 

Program  (EPWINDP) 

The  methodology  used  in  testing  EPWINDP  predicted  values 
against  observed  values  is  identical  to  that  used  in  the 
Atlantic  by  Jarrell^  (1981).  An  array  of  30  points  in  the 
eastern  North  Pacific  was  selected  (figure  1).  EPWINDP  values 
for  30  and  50  knots  were  calculated  at  12  hour  intervals  from 
the  effective  synoptic  time  of  the  Naval  Western  Oceanography 
Center  (NWOC)  forecasts  for  the  1980  season.  Since  most  of 
these  30  points  are  not  observing  stations,  actual  verifying 
winds  were  not  generally  available.  Consequently  a  verifying 
"warning  time"  probability  greater  than  50%  constituted  a 
verifying  strike. 

Tables  1,  2,  3  and  4  compare  the  expected  to  the  observed 
occurrences  of  30  and  50  knot  winds.  Predictions  are  associ¬ 
ated  with  percentage  groups  of  increasing  width,  <  1/2%,  1/2 
to  1  1/2%,  .  .  .  etc.  Time  integrated  probabilities  were 
verified  only  if  a  continuous  record  was  available  over  the 
entire  time  period. 
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gure  1.  Selected  test  points  (30)  in  the  Eastern  Pacific 
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Table  1.  Instantaneous  Probabilities. 

50  kt  winds  -  Expected  versus  Observed. 


30  KT 

24  Hr  Hr  72  Hr 


A<P<B 

N 

E 

0 

N 

E 

0 

N 

E 

0 

<  4% 

4978 

0 

5** 

3604 

1 

4** 

2296 

1 

1 

4-  14 

106 

1 

5** 

203 

2 

4 

258 

2 

2 

14  -  34 

75 

2 

4 

176 

4 

10* 

221 

5 

12 

34-  74 

80 

4 

8 

137 

7 

13* 

147 

7 

8 

74  - 154 

70 

8 

12 

93 

10 

11 

67 

7 

7 

154  -  314 

56 

12 

16 

47 

10 

9 

11 

3 

2 

314  -  634 

31 

13 

14 

0 

0 

0 

0 

0 

0 

>  634 

4 

3 

2 

0 

0 

0 

0 

0 

0 

ALL 

5400 

43 

66** 

4260 

34 

51* 

3000 

25 

32 

Table  2.  Instantaneous  Probabilities. 

30  kt  winds  -  Expected  versus  Observed 


*Slgnlf leant  at  5Z  level. 

**SlKniflcant  at  5%  level  even  after  ad.lustment  for  interdependence 
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Table  3. 

Time  Integrated  Probabilities. 

50  kt  winds  -  Expected  versus  Observed 
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2 

143 

1 

0 

li  -  3i 

63 

1 

2 

111 

3 

4 

103 

2 

1 

3i-  1\ 

79 

4 

4 

125 

7 

8 

134 

7 

5 

7i-15i 

93 

10 

10 

138 

15 

15 

164 

18 

18 

15i-  3li 

109 

24 

22 

145 

32 

38 

141 

30 

43 

31i  -  63j 

106 

48 

54 

99 

43 

51 

80 

35 

39 

>  63i 

34 

26 

26 

23 

17 

17 

16 

11 

11 

ALL 

5010 

114 

124 

3510 

118 

137 

2190 

104 

118 

Table  4.  Time  Integrated  Probabilities. 

30  kt  winds  -  Expected  versus  Observed. 

*Slgnlf leant  at  5%  level. 

^^Significant  at  52  level  even  after  adjustment  for  interdependence. 
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Significance  of  the  differences  between  the  expected  and 
the  observed,  as  discussed  in  previous  reports  (Jarrell^^, 
1981),  is  difficult  to  assess,  but  using  a  "t"  test,  agreement 
appears  to  be  very  good.  The  instantaneous  probabilities  show 
excellent  correlation  at  all  forecast  lengths  with  minor  under¬ 
forecasting  at  the  shorter  forecast  times  (24  and  48  hr).  The 
time  integrated  probabilities  displayed  a  slight  tendency 
toward  overforecasting  for  both  50  let  and  30  kt  winds.  The 
difference  between  the  expected  and  observed  occurrences  were 
rarely  at  a  significant  level  (see  note  accompanying  tables) 
and  the  overall  results  are  considered  statistically  sound. 

4 . 0  Operational  Products 

The  eastern  Pacific  wind  probability  program  will  be 
available  for  the  same  preselected  points  as  the  eastern 
Pacific  stri)ce  program.*  (The  EPSTRKP  program  is  actually  in¬ 
tegrated  into  the  EPWINDP  program. )  Probabilities  will  be 
given  in  two  modes,  instantaneous  and  time  integrated,  and  at 
0,  12,  24,  36,  48,  60  and  72  hours  after  the  warning  time.  The 
instantaneous  probabiity  will  be  the  probability  at  the  stated 
time  (i.e.,  12  hr)  and  the  time  integrated  probability  will  be 
summed  for  the  0  to  X  hour  time  interval  for  an  estimate  of  the 
probability  that  the  event  will  be  observed  within  that  period 
of  time. 


*Locations:  Acapulco,  Mazatlan,  Puerto  Vallarta,  La  Paz,  San 
Diego,  Hilo,  Honolulu,  Johnston  Is.,  Midway  Is. 


i 
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The  greatest  source  of  probable  error  for  the  EPWINDP 
program  will  be  erroneous  input  data.  An  internal  check  for 
unusual  motion  (expected  to  occur  only  5%  of  the  time  in 
nature)  will  be  made  and  suspect  motion  flagged.  The  user 
should  then  recheck  input  data  for  accuracy. 

When  the  forecast  track  approaches  a  land  mass,  the  fore¬ 
caster  should  be  aware  of  progrcun  bias.  This  should  be  minor 
for  seaward  approach  to  low  coastal  areas  or  over  smaller 
islands.  However,  in  other  cases  land  influences  will  appear 
as  rapid  decreases  in  the  instantaneous  wind  probabilities 
(especially  50  kt  winds)  near  forecast  landfall  time.  This 
will  bias  probabilities  -.overstate  them  for  inland  sites  and 
understate  them  for  coastal  sites.  Time  integrated  probabili¬ 
ties  will  be  less  biased.  This  problem  is  caused  by  wind  fore¬ 
casts  being  influenced  by  track  forecasts  where  landfall  is 
concerned.  A  bad  track  forecast  may  cause  a  bad  wind  forecast. 
This  was  not  accounted  for  either  in  development  nor  testing; 
hence  the  test  results  simulate  expected  actual  operational 
results  and  some  of  the  minor  disparities  between  expected  and 
observed  occurrences  no  doubt  stem  from  this. 
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5.0  Sxunmar 


The  wind  probability  model  for  the  eastern  North  Pacific 
is  largely  based  on  the  strike  probability  progreun  for  the 
eastern  North  Pacific.  Features  of  wind  probability  programs 
for  the  Atlantic  and  western  North  Pacific  have  been  incorpo¬ 
rated.  All  of  these  programs  are  operational  or  soon  to  be 
operational.  There  are  only  minor  modifications  to  those  com¬ 
ponents  adapted  herein.  Test  results  of  the  eastern  North 
Pacific  wind  probability  program  (EPWINDP)  demonstrated  excel¬ 
lent  agreement  between  expected  and  observed  results. 
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